function [v1_new, v2_new] = updateVectors(f, a, x, phi, n, v1, v2)
    % updateVectors 利用有限差分法计算函数 f 在 (x, phi) 处的雅可比矩阵
    % f: 函数句柄, 调用形式为 [x_next, phi_next] = f(a, x, phi, n)
    % a, x, phi, n: 当前参数和状态
    % v1, v2: 当前扰动向量
    %
    % 输出:
    %   v1_new, v2_new: 更新后的扰动向量
    
    h = 1e-8;  % 小扰动步长
    
    % 计算在 (x, phi) 处的 f 值
    [fx, fphi] = f(a, x, phi, n);
    F = [fx; fphi];
    
    % 对 x 方向做扰动
    [fx_x, fphi_x] = f(a, x + h, phi, n);
    F_x = [fx_x; fphi_x];
    dfdx = (F_x - F) / h;
    
    % 对 phi 方向做扰动
    [fx_phi, fphi_phi] = f(a, x, phi + h, n);
    F_phi = [fx_phi; fphi_phi];
    dfdphi = (F_phi - F) / h;
    
    % 构造雅可比矩阵
    J = [dfdx, dfdphi];
    
    % 更新扰动向量
    v1_new = J * v1;
    v2_new = J * v2;
end
